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Abstract 

We study transformations of 2-, 4- and 6-bit numbers in interactions 
between traveling and stationary localizations in the Spiral Rule reaction- 
diffusion cellular automaton. The Spiral Rule automaton is a hexagonal 
ternary-state two-dimensional cellular automaton - a finite-state machine 
imitation of an activator-inhibitor reaction-diffusion system. The activa- 
tor is self-inhibited in certain concentrations. The inhibitor dissociates 
in the absence of the activator. The Spiral Rule cellular automaton has 
rich spatio-temporal dynamics of traveling (glider) and stationary (eater) 
patterns. When a glider brushes an eater the eater may slightly change 
its configuration, which is updated once more every next hit. We encode 
binary strings in the states of eaters and sequences of gliders. We study 
what types of binary compositions of binary strings are implementable 
by sequences of gliders brushing an eater. The models developed will be 
used in future laboratory designs of reaction-diffusion chemical computers. 



Keywords: cellular automata, reaction-diffusion computing, gliders, collision- 
based computing 

1 Introduction 

Many physical, chemical and biological spatially extended non-linear systems 
exhibit a wide range of stationary and mobile localizations: solitons, kinks, 
breathers, excitons, defects and wave-fragments. The localizations can be used 
to transmit and transform information, and ultimately to perform computa- 
tion PQ. A unit of information, such as the value of a Boolean variable, is 
decoded into presence (logical truth) or absence (logical false) of a localiza- 
tion in some specified site of space at a specified moment of time. When 
two localizations (representing the values of two logical variables) collide, they 
change their trajectories (or annihilate, reproduce, or change their shape). The 
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new trajectories of the localizations encode the values of some logical function 
over the two variables. This is how most collision-based computing devices 

work is uni nai m m na eui si]- 

The collision-based, or free-space, computing devices typically do not have 
wires and — in principle - are not supposed to use any other stationary com- 
ponents to perform computation. Any point of the computing media can act 
as a wire, a trajectory of a traveling localization can be seen as a momentary 
wire. Any site where two or more localizations collide is a logical gate. Thus 
space can be used efficiently and nothing is wasted. However, there is a price 
to pay. Initial positions and launch time of the traveling localizations should be 
precisely specified: one wrong time step destroys the whole computing scheme. 
The need for perfect timing is the weakest point of collision-based computing 
and the subject of constant criticism by "classical" computation schools. 

Can we do without perfect timing? Asynchronous cellular-automaton based 
computers do pretty well I15[ 112] by using predetermined wires and valves. 
In the present paper we are trying to combine pure collision-based computing 
ideas (gliders only) with stationary architectures (breather-like localization) to 
implement computing schemes with relaxed timing. We choose the reaction- 
diffusion Sprial Rule cellular automaton [18j E 13 HI] as a testbed for ideas of 
asynchronous collision-based computing. 

What is our rationale behind selecting the Spiral Rule automaton to study 
novel concepts of the collision-based computing? The Spiral Rule cellular au- 
tomaton |18j plays a unique role in unconventional computing. On the one hand, 
this is a simple ternary state hexagonal automaton with Conway's Game of Life 
type of behavior: it has gliders, still lives and eaters, and glider. Therefore it is 
very suitable for experimenting with collision-based computing schemes. On the 
other hand, the Spiral Rule automaton is a unique discrete model of a non-linear 
reaction-diffusion chemical system with an activator and inhibitor. The gliders 
and glider guns in the Spiral Rule automaton are analogues of excitation wave- 
fragments and generators of wave-fragments in a light-sensitive sub-excitable 
Belousov-Zhabotinsky medium [9] . This means that prototypes of computing 
schemes designed in the Spiral Rule automaton can then be almost straight- 
forwardly implemented in chemical laboratory prototypes of reaction-diffusion 
computers. 

The paper is structured as follows. The Spiral Rule reaction-diffusion cellular 
automaton is defined in Sect.[2j Encoding of two- and four-bit binary strings in 
states of an eater and transformation of the strings by gliders are presented in 
Sect. [3] Section [4] studies compositions of six-bit binary strings and speculates 
on the algebraic properties that the compositions impose on a set of six-bit 
digital numbers. We envisage future developments in Sect. [5] 

2 Spiral rule reaction-diffusion cellular automa- 
ton 

A detailed description, including the background and rationale, of the Spiral 
Rule reaction-diffusion cellular automaton are provided in [18j|4j|3j[l6], therefore 
here we just give a self-consistent outline. The Spiral Rule automaton is a 
discrete approximation of an abstract reaction-diffusion chemical system. The 
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system has three reactants — activator A, inhibitor / and substrate S — and 
the system is governed by the following set of reactions [TH1 H]: 
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The system is non- linear. The activator is suppressed by the inhibitor when 
there are small concentrations of the inhibitor and threshold concentrations of 
the activator. When the inhibitor has a critical concentration (value 2) the 
inhibitor and activator dissociate producing the substrate. 

In the cellular-automaton model each rcactant is represented as a cell state. 
Each cell i of a two-dimensional hexagonal lattice updates its state by the 
following rule: 

where a p (x) t is the number of cell x's neighbors with cell-state p £ {I, A, S} 
at time step t. Cells update their states synchronously in discrete time-steps. 
The cell-state transition rule can be compactly represented by matrix P = (pij), 
where < i < j < 7, < % + j < 7, p tj £ {I, A,S} [3J [T5] as follows: 

P(T2 (x) t (Ti(x) t • 

The matrix has the following structure: 
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The cell-state transition matrix could be simplified further because the in- 
hibitor states underlined above are 'wildcard' entries [TS]. In [J] we demon- 
strated that each entry of the matrix P corresponds very well with phenomena 
in reaction-diffusion systems: diffusion of the activator, suppression of the acti- 
vator by the inhibitor, self-inhibition of the activator etc. 

Any initially random configuration of the Spiral Rule automaton evolves 
to a dissipative configuration comprised of spiral glider guns (Fig. [TJ, mobile 
localizations (gliders) and stationary localized clusters of activator and inhibitor 
states (eaters). The principle, and of minimal size, localizations are shown in 
Fig. [2] A glider always has one activator state tailed by a few inhibitor states 
(Fig. [21 a-c), pretty much as excitation wave-fragment in sub-excitable chemical 
medium 2 . The eater (Fig.[2ji) — in its minimal symmetric form [J] — consists 
of six activator states surrounded by seven inhibitor states (including the central 
state of the eater) . 
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Figure 1: The full cycle of the high-frequency spiral glider gun. The "Core" of 
the gun is shown rotating (similarly to spiral waves in a sub-excitable Belousov- 
Zhabotinsky) medium [2]), but can rotate in either direction. The gun emits six 
streams of gliders. An alternative, low-frequency, spiral glider gun also emerges 
in the Spiral Rule QUE!]. 

4 




(a) t + 2 (b) t + 1 (c) t 




(d) 



Figure 2: Basic localizations in Spiral Rule. Mobile localizations, gliders types 
gl-g5, and the stationary localization, eater (d). Gliders are shown moving from 
right to left - time steps (c) to (a). Gliders gl, g2 and g5 change shape between 
alternate time-steps, gliders g4 and g5 arc asymmetric with both chiral forms 
present. Cell-state / (inhibitor) is shown by a circle, cell-state A (activator) by 
a solid disk, and cell-state S (substrate) by a dot. 
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The Spiral Rule supports a rich population of other structures — compound 
gliders, mobile glider guns, compound eaters [TT1 [TH], in addition to the basic 
ones discussed. There is also another basic stable localization SL1 [17l[T8], which 
can modify gliders brushing past, but does not have "memory" in its outer shell 
like the eater discussed in the present paper. 

3 Transformation of two- and four-bit strings 

In [4] we discovered that substrate-sites between inhibitor-sites of the eater can 
be switched to an inhibitor-state by a colliding glider, or even a glider just 
brushing past. Thus the eater can play a role of a non- volatile memory device 
and also implement an asynchronous XOR gate [I]. What if we consider the 
states of several substrate-sites at once? Then gliders brushing the eater can 
change these states recursively and thus modify the string content represented 
by the states. 

Let us encode two substrate-states in the northern part of the eater as X\ and 
X2 (Fig. |3^i). If the cell Xi takes the inhibitor state it represents a bit-up state 
Xi of the string (0:1X2), if the cell x, is in the substrate-state then it represents 
a bit-down state of the string (X1X2). We have four possible combination of 
cell-states xi and X2: SS, SI, IS, II, which corresponds to binary string 00, 
01, 10 and 11, respectively (Fig. [3]d). 

An example of a binary string transformation is shown in Fig. [4j Initially, 
the eater represents string 00, cells Xi and X2 are in the substrate state (Fig. [4^). 
We launch glider g\ traveling West (Fig. |1|d). When the glider brushes along 
the outer edge of the eater (Fig. |4^,)f-k) it causes a minor perturbation of the 
marginal cell-states. The result of the perturbation is that one substrate-cell Xi 
becomes an inhibitor (Fig. |4^,)k). The glider changes its state as well — g\ — > g$ 
— because of its interaction with the eater (Fig. |3Ji)j-o). 

All configurations, including those modified by gliders, are still configura- 
tions, they arc not evolving by themselves but can only be changed again by 
another brushing glider. When we apply a sequence of transformation we do not 
have to worry about the precise timing of glider generation or a time interval 
between gliders (as far as the gliderphase is preserved). Therefore, eater-glider 
based transformations can be classified as asynchronous. 

Finding 1 The following transformations of two-bit strings, encoded below in 
digital numbers, are implemented by gliders and eaters in the Spiral Rule cellular 
automaton: 

L\ : L 2 : 

0^0 ^ 2 ^ L 3 :0^1^2^3^0 

1 -> 3 -> 1 1—^1 L 4 :0^3^2^1^0 
2^2 3^3 



Assuming states '0' and T' are Boolean values, transformations L\ ■ ■ ■ L3 can 
be expressed in the following manner: 

Li(xi,Xi_i) = (x.i V yi A (x.i © Xi_i),Xi_i V yi A x^_i) (1) 
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Figure 3: Encoding two-bit strings in the states of eaters. The state of cells x\ 
and £4 for a single cater encoding (a) and cells x\ . . . X4 for the eater couple (c) 
encoding represent bits of strings (212:2) and (X1X2X3X4) respectively. Cell-state 
/ (inhibitor) is shown by a circle, cell-state A (activator) by a solid disk, and 
cell-state 5* (substrate) by a dot. Cell Xi in state I represents T', cell in state S 
represents '0'. Examples of encoding strings 00, 01, 10, 11 in a single eater are 
shown in (b). Examples of encoding strings 1110 and 0110 in the eater couple 
are shown in (d). 
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Figure 4: Transformation 00 — > 10 implemented by glider g\. 
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L 2 (x i ,x i - 1 ) = (xi V yi A (xi © Xi-x), Xi-i V yi A (2) 



i 3 (x i ,x i _i) = (x.i V j/i A (x.i © £j_i),Xi_i V y« A Xj-i) (3) 

L 4 (x 4 ,x 4 _i) = (if V t/i A (xj © Xj_i),Xi_i V A x~T). (4) 

Finding 2 One can implement increment modulo 4 using transformation L 3 
and decrement modulo 4 using transformation L% 

This follows directly from Finding [I] and Boolean representation ^ . 

Similarly by positioning two eaters close to each other — the distance between 
central states of eaters is 4 — we can change four substrate-states just with one 
glider. The positions of substrate-states x% . . . X4 used for encoding are shown 
in Fig. [3p; examples of encoding four-bit binary sequences are demonstrated in 
Fig.p. 

Example of transformation 0000 — > 1100 is shown in Fig. [5] and transforma- 
tion 1100 ->■ 0010 in Fig. [6] 

Finding 3 The following transformations of four-bit strings, encoded below in 
digital numbers, are implemented by gliders and eaters in the Spiral Rule cellular 
automaton: 
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The exact conditions C(Tj) or C(Zj) of transformations Tj and 1^ are deter- 
mined by two parameters (g,d), where a type g glider brushes the eater, and 
the initial distance d between the eater and the glider is even. 

Finding 4 The conditions of the transformations Li and Ti, i = 1 . . .4 are as 
follows: 

• C(Li) — C(Xi) = {g 1 ,d'), g' =53,34 and d'=even, odd. 



9 




(d) (c) (f) 




(j) (k) (1) 




(P) 



Figure 5: Transformation 0000 — > 1100 implemented by glider g 3 . 
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• C{L2) — C(T2) — (g',df), g' — g\ and dl —even, odd. 

• C(L 3 ) = C(T 3 ) = {(g 2 , even), (g 7 , even), (g 5 , odd)} 

• C(L 4 ) = C(T 4 ) = {(.92, odd), (g 7 , odd), (g 5 , even)}. 

Finding 5 Transformations T\ and Ti split the set {0, • • • , 15} into stable sub- 
sets of exclusively even numbers and subsets of exclusively odd numbers. 

This is an outcome of the experimental implementation of the transforma- 
tions Ti and T 2 , and analysis of chains generated by the transformations in 
Finding [3] 

Finding 6 By using a series of transformations T\, . . . , T 4 one can transform a 
number from {0, . . . , 15} to any other number from {0, . . . , 15} in at most three 
steps. 

To prove this statement we compute a distance matrix D = (drj)o<z,j<i5, 
where dij represents the minimal number of transformations T\ , . . . T4 that one 
needs to apply to get from number i to number j. The distance matrix has the 
following form: 

/3221231221321223 v 
/ 32212 3 1221321223 \ 

32132 3 3 11 3 322123 

3213233113322123 

1332312332132331 

1332312332132331 

2132122332212312 
n _ 21 3 21223 3 2212312 
^ — 2132122332212312 

2132122332212312 

1332312332132331 

1332312332132 3 3 1 

32132 3 311332212 3 

32132 3 3113321223 
\3221231221321223 / 
^3221231221321223 / 

For example, to get 8 from 7 we must apply T\ twice and then T3 once: 
Ts(T\(T\(J))) = 8. From D we can also tell that one needs 34 to generate a 
glider, which brushes the eater to generate a monotonous ordered sequence of 
numbers (0,1,..., 15). 

4 Six bit coding 

Let us now consider the composition xoy of two six bit binary strings x = 
(xi-xg) and y = (jjvyo). The string x is represented by states of six cells x\ ■ ■ ■ xg 
surrounding the eater, as shown in Fig. [7^. The string y is represented by gliders 
brushing the eater (as discussed in the previous sections). The trajectories of 
gliders are indicated in Fig. [7] yi = 1 if there is a glider traveling along trajectory 
yi (Fig.Qa. Examples of encoding binary strings in the eater's states are shown 
in (Fig.|7|b. 

Gliders representing string y are sent in the order y± , • • • , y@ . Values of 
Boolean variables yi, i — 1, • • • ,6, are represented by the presence or absence 
of a glider g, traveling along trajectory yi. A minimum interval between two 
subsequent gliders is selected so the next glider is sent after the current glider has 
already passed the eater. There is no maximal interval between the gliders. The 
exact structure of composition xoy is determined by transformations L\ ■ ■ ■ L4. 



12 



yi 




(a) (b) (c) 



Figure 7: Encoding six-bit strings in states of the eaters, (a) States of cells 
X\ ■ ■ -Xq represent bits of strings (xi • • • Xq). Cell-state / (inhibitor) is shown 
by a black disk, cell-state A (activator) by a circle, and cell-state S (substrate) 
by a dot. Cell x, in state L represents '1', cell in state S represents '0'. (b)- 
(c) Examples of encoding strings 101100 (a) and 001101 (b). Vectors y% ■ ■ ■ ye 
in (a) show which direction gliders representing input string y± ■ ■ ■ ye travel. 



For example, if the composition xoy is based on L4 then the result xVxoy 
can be calculated as follows: 

x[ = xi V yi A (xi © x 6 ) 
x'e = x 6 V yi A x~e~ 



and for i = 2, . . . , 6 



x- = Xi V yi A (xl © Xi_i) 
x- = Xi_a V yi A x^TI 



We represent the composition xoy implemented by the eater and glider g z , 
z = 1, • • • , 6, by matrix M{Li) = (?riij)o<ij<63- Each entry niij of the matrix 
is an integer from {0, • • • ,63}. The entry rriij is a digital representation of the 
string x<-xoy, and its indexes i and j are digital representations of strings x 
and y, respectively. 

Matrices of the compositions xoy are shown in Appendix (Fig. [8 11 1. Struc- 
tures of the matrices do not indicate any immediate implications for practical 
applications. There is however a range of properties which may be useful in 
further implementations of arithmetical operations with gliders and eaters. 

Finding 7 Let k(L{), i = 1, ... ,4, be a ratio of commutative pairs in xoy 
determined by transformation Li, and a(Li) a ratio of associative triples, then 
the following hierarchy of compositions can be drawn: 

k(L 3 ) > k(Li) > k(L 4 ) > k(L 2 ) 

a(L 2 ) > a(L 3 ) > a(Li) > a(L±) 

This follows from straightforward analysis of the matrices M(Li), . . . , M(L^) 
in Tab. CD 

Finding 8 For compositions determined by transformations L\ and L 2 we have 
x o x = x for any x from {0, . . . , 63}. 
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Li 



n{Lj) a{L % ) 



Li 
L 2 
L, 
Li 



0.33 0.79 

0.11 0.85 

0.35 0.81 

0.15 0.63 



Tabic 1: Ratios of commutative pairs «;(•) and associative triples a(-) in 



The following finding relates to identities and absorbing elements (nulls) of 
the compositions. 

Finding 9 Number 63 is the absorbing element of compositions determined by 
£3 and L4. The composition determined by Li has two left absorbing elements: 
and 63. The composition determined by L2 has eighteen left absorbing elements: 



21, 23, 27, 29, 31, 42, 43, 45, 46, 47, 53, 54, 55, 58, 59, 61, 62, 63. None of 



the composition has left or right identities. 

The structure of the first two rows of M{L\) shows that stationary and 
traveling localizations in the Spiral Rule automaton can implement at least one 
sensible arithmetical operation: transformation of arbitrary stream of numbers 
to streams of even or odd numbers. 

Finding 10 Let number a € {0, 1}, encoded in a state of the eater, and number 
b G {0, . . . , 31}, encoded by a stream of gliders traveling along trajectories as as 
per Fig. [?] then the eater-gliders system executes operation 2b + a; the result of 
the operation is encoded in state of the eater. 

The resultant state of the eater can be read and erased by additional stream 
of gliders, as have been already demonstrated in |3]. 

5 Discussion 

Using a cellular-automaton model [I] of a reaction-diffusion system we uncovered 
ways of implementing asynchronous operations over sets of binary strings. We 
employ a sequence of mobile localizations (gliders, excitation wave-fragments) 
to represent one binary string, and a stationary localization (eater, breather, 
standing wave) to represent another string in the composition. The operations 
are asynchronous in the sense that there is no precise time for a glider to collide 
with the eater, and there is no fixed time interval between gliders in a glider 
train representing a binary string. The order in which gliders hit certain parts 
of the eater is nevertheless important and must be executed precisely during 
the computation. 

We found that by colliding gliders with eaters we can implement transfor- 
mations of sets of two-, four- and six-bit integers, which demonstrates collision 
based computation in the Spiral Rule. The methods could find applications for 
the design arithmetical chips. 

We did not discover any sophisticated operations on digital numbers. Most 
impressive so far were increments and decrements modulo four and subdivision 
of the digital number set onto prime and even numbers. We found operations 



M{L\), . . . , M(L 4 ). 
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on a set of six-bit integers which are rich yet irregular and require detailed 
study. Nevertheless are we convinced that by combining operations (based on 
L\, . . . ,Li) one can implement typical arithmetical operations. This will be 
subject of further analysis. 

What are the potential practical outcomes of our designs and findings? As 
we demonstrated in [14] . the Spiral Rule automaton is an ideal discrete ap- 
proximation of excitation dynamics in a light-sensitive sub-excitable chemical 
medium on a patterned environment. Gliders are discrete analogs of localized 
wave-fragments. This does not mean that our computing schemes can be imple- 
mented directly in an excitable chemical medium, not least because there are 
no excitation-wave analogues of eaters. More feasible implementations would 
be in other substrates, possibly at the nano-scales as e.g. excitation patterns 
in networks of single-electron oscillators jS], excitons in mono-molecular layers 
and breathers in polymers and crystals pQ, and dynamically evolving molecular 
circuits [6]. 
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